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One explanation for the impressive recent boom in network theory might be that it provides a promising tool for an 
understanding of complex systems. Network theory is mainly focusing on discrete large-scale topological structures 
rather than on microscopic details of interactions of its elements. This viewpoint allows to naturally treat collective 
phenomena which are often an integral part of complex systems, such as biological or socio-economical phenomena. 
Much of the attraction of network theory arises from the discovery that many networks, natural or man-made, seem 
to exhibit some sort of universality, meaning that most of them belong to one of three classes: random, scale-free 
and small-world networks. Maybe most important however for the physics community is, that due to its conceptually 
intuitive nature, network theory seems to be within reach of a full and coherent understanding from first principles. 

Networks are discrete objects made up of a set of nodes which are joint by a set of links. If a given set of N nodes 
is linked by a fixed number of links in a completely random manner, the result is a so-called random network, whose 
characteristics can be rather easily understood. One of the simplest measures describing a network in statistical terms 
is its degree distribution, p{k), (see box 1). In the case of random networks the degree distribution is a Poissonian, 
i.e., the probability (density) that a randomly chosen node has degree k is given by p{k) = — , where A = is 
the average degree of all nodes in the network. However, as soon as more complicated rules for wiring or growing of 
a network are considered, the seemingly simple concept of a network can become rather involved. In particular, in 
many cases the degree distribution becomes a power-law, bare of any characteristic scale, which raises associations to 
critical phenomena and scaling phenomena in complex systems. This is the reason why these types of networks are 
often called complex networks. A further intriguing aspect of dynamical complex networks is that they can potentially 
provide some sort of toy-model for 'nonergodic' systems, in the sense that not all possible states (configurations) 
are equally probable or homogeneously populated, and thus can violate a key assumption for systems described by 
classical statistical mechanics. 

Over the past two decades the concept of nonextensive statistical mechanics has been extremely successful in 
addressing critical phenomena, complex systems and nonergodic systems 0, Q . Nonextensive statistical mechanics is 
a generalization of Boltzmann-Gibbs statistical mechanics, where entropy is defined as 

1 - f°° dk \v(k)]'^ 

Sq = with the g = 1 limit Si ^ Sbg = - dkp{k) \np{k) (1) 

9^1 Jo 

where BG stands for Boltzmann-Gibbs. If - in the philosophy of the maximum entropy principle - one extremizes Sq 
under certain constraints, the corresponding distribution is the ^-exponential (see box in editorial paper by C. Tsallis 
and J. P. Boon). Another sign for the importance and ubiquity of (/-exponentials in nature might be due to the fact 
that the most general Boltzmanfactor for canonical ensembles (extensive) is the g-exponential, as was recently proved 
in Pi . Given the above characteristics of networks and the fact that a vast number of real- world and model networks 
show asymptotic power-law degree distributions, it seems almost obvious to look for a connection between networks 
and nonextensive statistical physics. 

Since the very beginning of the recent modeling efforts of complex networks it has been noticed that degree distri- 
butions asymptotically follow power-laws 0, or even exactly g-exponentials 5j. The model in |^] describes growing 
networks with a so-called preferential attachment rule, meaning that any new node being added to the system links 
itself to an already existing node i in the network with a probability that is proportional to the degree of node i, i.e. 
piink fx ki. In this model was extended to also allow for preferential rewiring. The analytical solution to the model 
has a g-exponential as a result, with the nonextensivity parameter q being fixed uniquely by the model parameters. 

However, it has been found that networks exhibiting degree distributions compatible with g-exponentials are not 
at all limited to growing and preferentially organizing networks. Degree distributions of real-world networks as well 
as of models of various kinds seem to exhibit a universality in this respect. In the remainder we will review a small 
portion of the variety of networks which potentially have a natural link to non-extensive statistics. Even though there 
exists no complete theory yet, there is substantial evidence for a deep connection of complex networks with the q ^ 1 
instance in nonextensive statistical mechanics. 

Recently in preferential attachment networks have been embedded in Euclidean space, where the attachment 
probability for a newly added node is not only proportional to the degrees of existing nodes, but also depends on 
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the Euclidean distance between nodes. The model is realized by setting the linking probability of a new node to an 
existing node i to be punk ^ ki/rf [a > 0), with being the distance between the new node and node i; a = 
corresponds to the model in Q which has no metrics. In a careful analysis the degree distributions of the resulting 
networks have been shown to be g-exponentials with a clear a-dependence of the nonextensivity parameter q. In the 
large a limit, q approaches one, i.e., random networks are recovered in the Boltzmann-Gibbs limit. 

In an effort to understand the evolution of socio-economic networks, in a model was proposed that builds upon 0| 
but introduces a rewiring scheme which depends on the internal network distance between two nodes, i.e., the number 
of steps needed to connect the two nodes. The emerging degree distributions have been subjected to a statistical 
analysis where the hypothesis of g-exponentials could not be rejected. 

A model for nongrowing networks which was recently put forward in |^ also unambiguously exhibit g-exponential 
degree distributions. This model was motivated by interpreting networks as a certain type of 'gas' where upon an 
(inelastic) collision of two nodes, links get transfered in analogy to the energy-momentum transfer in real gases. In this 
model a fixed number of nodes in an (undirected) network can 'merge', i.e., two nodes fuse into one single node, which 
keeps the union of links of the two original nodes; the link connecting the two nodes before the merger is removed. 
At the same time a new node is introduced to the system and is linked randomly to any of the existing nodes in the 
network . Due to the nature of this model the number of links is not strictly conserved - which can be thought of 
as jumps between discrete states in some 'phase space'. The model has been further generalized to exhibit a distance 
dependence as in 0, however not being Euclidean but internal distance. Again, the resulting degree distributions 
have q-exponential form. In Fig. 1 we show a snapshot of this type of network pars pro toto for the many models 
exhibiting g-exponential degree distributions. The corresponding (cumulative) degree distribution is shown in Fig. 
2 in log-log scale, clearly exhibiting a power-law. Figure 3 shows g-logarithms of the degree distribution for several 
values of q. It is clear from the correlation coefficient of the q-logarithm with straight lines (inset) that there exists 
an optimal value of q, which makes the g-logarithm a linear function in fc, showing that the degree distribution is a 
^-exponential. _ 

A quite different approach was taken in 10] where an ensemble interpretation of random networks has been adopted, 
motivated by superstatistics (lli] . Here it was assumed that the average connectivity A in random networks is fluctu- 
ating according to a distribution n(A), which is sometimes associated with a 'hidden- variable' distribution. In this 
sense a network with any degree distribution can be seen as a 'superposition' of random networks with the degree 
distribution given by p{k) = dAn(A) ^ ■ In 0| it was shown as an exact example, that a power-law func- 
tional form of n(A) leads to degree distributions of Zipf-Mandelbrot form, p{k) oc (fep^^-)-, , which is equivalent to a 
g-exponential with an argument of k/n and given the substitutions, k = (1 — q)kQ and q = I + I/7. Very recently 
a possible connection between small-world networks and the maximum Tsallis-entropy principle, as well as to the 
hidden variable method 10], has been noticed in 12]. 

In yet another view of networks from a physicist's perspective, networks have recently been treated as statistical 
systems on a Hamiltonian basis. It has been shown that these systems show a phase transition like behavior |13], along 
which networks structure changes. In the low temperature phase one finds networks of 'star' type, meaning that a few 
nodes are extremely well connected resulting even in a discontinuous p{k); in the high temperature phase one finds 
random networks. Surprisingly, for a special type of Hamiltonians networks with g-exponential degree distributions 
emerge right in the vicinity of the transition point [l^ . 

While a full theory of how complex networks are connected to q ^ 1 statistical mechanics is still missing, it 
is almost clear that such a relation should exist. It would not be surprising if an understanding of this relation 
would arise from the very nature of networks, being discrete objects. More specifically, it has been conjectured for 
nonextensive systems that the microscopic dynamics does not fill or cover the space of states (e£. F-space (6A^ 
dimensional phase space) for Hamiltonian systems) in a homogeneous and equi-probable manner |3]. This possibly 
makes phase space for nonextensive systems look like a network itself, in the sense that in a network not all possible 
positions in space can be taken, but that microscopic dynamics is restricted onto nodes and links. In this view the 
basis of nonextensive systems could be connected to a network like structure of their 'phase space'. It would be 
fantastic if further understanding of network theory could propel a deep understanding of nonextensive statistical 
physics, and vice versa, making them co-evolving theories. 

BOX 1: Some network measures 

The degree ki of a particular node i of the network is the number of links associated with it. If links are directed 
they either emerge or end at a node, yielding the diction of out- or in-degree, respectively. The degree distribution 
p{k) is the probability for finding a node with degree k in the network. In (unweighted) networks the degree 
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distribution is discrete and often reads, p{k) = p{l)eq with k > being some characteristic number of hnks. Apart 
from the degree distribution, important measures to characterize network topology are the clustering coefficient Ci, 
and the neighbor connectivity fc„n- The clustering coefficient measures the probability that two neighboring nodes 
of a node i axe also neighbors of each other, and is thus a measure of cliquishness within networks. The neighbor 
connectivity is the average degree of all the nearest neighbors of node i. When plotted as a function of k, a nontrivial 
distribution of the average of c allows statements about hierarchic structures within the network, while fc„„ serves as 
a measure of assortativity. 



[1] C. Tsallis. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys., 52:479-487, 1988. 

[2] M. Gell-Mann and C. Tsallis. Nonextensive Entropy - Interdisciplinary Applications. (Oxford University Press, New York, 
2004); C. Tsallis, M. Gell-Mann, and Y. Sato. Asymptotically scale-invariant occupancy of phase space makes the entropy 
Sq extensive. Proc. Natl. Acad. Sc. USA, in press, 2005; cond-mat/0502274 

[3] R. Hanel and S. Thurner. Derivation of power-law distributions within standard statistical mechanics. Physica A, 351:260- 
268, 2005. 

[4] A.-L. Barabasi and R. Albert. Emergence of scaling in random networks. Science, 286:509- 512, 1999. 
[5] R. Albert and A.-L. Barabasi. Topology of evolving networks: local events and universality. Phys. Rev. Lett., 85:5234-5237, 
2000. 

[6] D. J. B. Soares, C. Tsallis, A. M. Mariz, and L. R. da Silva. Preferential attachment growth model and nonextensive 

statistical mechanics Europhys. Lett., 70:70-76, 2005. 
[7] D. R. White, N. Kejzar, C. Tsallis, J. D. Farmer, and S. White. A generative model for feedback networks. 

cond-mat/0508028 

[8] S. Thurner and C. Tsallis. Nonextensive aspects of self-organized scale-free gas-like networks. Europhys. Lett., 72:197-203, 
2005. 

[9] B. J. Kim, A. Trusina, P. Minnhagen, and K. Sneppen. Self organized scale-free networks from merging and regeneration. 

Eur. Phy J. B, 43:369-372, 2005. 
[10] S. Abe and S. Thurner. Complex networks emerging from fluctuating random graphs: analytical formula for the hidden 

variable distribution. Phys. Rev. E, 72:036102, 2005. 

[11] C. Beck and E. G. D. Cohen. Superstatistics. Physica A, 322:267-275, 2003. 

[12] H. Hasegawa. Nonextensive aspects of small-world networks. Physica A, in press, 2005; |cond-mat/0506301 | 

[13] G. Palla, L Derenyi, L Farkas, and T. Vicsek. Statistical mechanics of topological phase transitions in networks. Phys. 

Rev. E, 69:046117, 2004. 

[14] C. Biely and S. Thurner. Statistical mechanics of scale-free networks at a critical point: complexity without irreversibility. 
|cond-mat/0507670, 



4 




FIG. 1; Snapshot of a non-growing dynamic network with g-exponential degree distribution for A'^ = 256 nodes and a linking 
rate of f = 1, for details see |^|^- The shown network is small to make connection patterns visible. 




FIG. 2: Log-log representation of (cumulative) degree distributions for the same type of network with various system sizes A'^ 
and f — 8. The distribution functions are from individual network configurations without averaging over identical realizations. 
On the right side the typical exponential finite size effect is seen. 
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FIG. 3: g-logarithm of the (cumulative) distribution function from the previous figure as a function of degree k. Clearly, there 
exists an optimum q which allows for an optimal linear fit. Inset: Linear correlation coefficient of In, P{> k) and straight lines 
for various values of q. The optimum value of q is obtained when In, P(> k) is optimally linear, i.e., where the correlation 
coefficient has a maximum. A linear In, means that the distribution function is a g-exponential; the slope of the linear function 
determines hi. In this example we get for the optimum q — 1.84, which corresponds to the slope 7 = 1.19 in the previous plot. 



